arXiv:l508.07138v2 [cond-mat.mes-hall] 4 Sep 2015 


Hedgehog spin texture and competing orders associated with strains on the surface of 

a topological crystalline insulator 

Cheng-Yi Huang, 1 Hsin Lin, 2, 3, [2] Yung Jui Wang, 4 Arun Bansil, 4 and Wei-Feng Tsai 1 ’^ 

department of Physics, National Sun Yat-sen University, Kaohsiung 80424, Taiwan 
2 Centre for Advanced 2D Materials and Graphene Research Centre, 

National University of Singapore, Singapore 117546 
department of Physics, National University of Singapore, Singapore 117542 
4 Department of Physics, Northeastern University, Boston, Massachusetts 2115, USA 

(Dated: September 7, 2015) 

We have investigated spin reorientation phenomena and interaction driven effects under the pres¬ 
ence of applied strains on the (001) surface of Pbi_ a ,Sn a: (Te, Se) topological crystalline insulators, 
which host multiple Dirac cones. Our analysis is based on a four-band k ■ p model, which captures 
the spin and orbital textures of the surface states at low energies around the X and Y points, 
including the Lifshitz transition. Even without breaking the time-reversal symmetry, we find that 
certain strains which break the mirror symmetry can induce hedgehog-like spin texture associated 
with gap formation at the Dirac points. The Chern number of the gapped surface ground state is 
shown to be tunable through the interplay of strains and a perpendicular Zeeman field. We also 
consider effects of strain in the presence of interactions in driving competing orders, and obtain the 
associated phase diagram at the mean-field level. Potential applications of our results for low power 
consuming electronics are discussed. 


I. INTRODUCTION 

A topological insulating (TI) phase is a new state 
of quantum matter featuring massless Dirac-like bound¬ 
ary states whose robustness is guaranteed by the time- 
reversal symmetry (TRS)2®I In this rapidly growing 
field, search for new topological phases has recently 
turned to extending the consideration of symmetry- 
protected states to include non-spatial symmetries 5 -*® 
such as, the particle-hole and chiral symmetries as well as 
the spatial symmetries' 11 . In particular, the nontrivial 
band topology can be shown to be protected by certain 
crystal symmetries, leading to the new class of TIs called 
topological crystalline insulators (TCis(P^b-4 Like the 
TIs, a three dimensional (3D) TCI is also predicted to 
host metallic surface states on surfaces which preserve 
appropriate crystal symmetrieif 15 ! 

To date, the semiconducting Pbi_ x Sn x (Te, Se) al¬ 
loys are the only experimentally realized TCI materials, 
whose topological nature has been verified through direct 
observation of Dirac-like surface states via angle-resolved 
photoemission spectroscopy (ARPES) experiment^ 1618 . 
These gapless surface states are protected by the mir¬ 
ror symmetry with respect to the (110) or (110) lat¬ 
tice plane. The characteristic features of the associated 
Dirac states such as linear dispersion, Lifshitz transi¬ 
tion, spin/orbital texture, among others, have been ex¬ 
amined via scanning tunneling microscop y/spe ctroscopy 
(STM/STS) and transport measurements 16 19 23 . More¬ 
over, it has been found that, below a critical temperature, 
which depends on chemical composition, the cubic lattice 
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structure can be distorted by strains, resulting in an or¬ 
thorhombic or a rhombohedral structure 24 . Remarkably, 
in a recent STM experiment by Okada et alW^, two un¬ 
expected Landau levels have been observed for the (001) 
surface of Pbi_ x Sn x Se under a perpendicular magnetic 
field. The presence of these two extra levels is believed 
to be associated with the gap opening of the two surface 
Dirac cones induced by a ferroelectric-like lattice distor¬ 
tion, which breaks the corresponding mirror symmetry. 
It is clear that it is important to understand effects of 
strains on the electronic structure of TCIs in order to 
gain a deeper handle on the nature and origin of their 
topological states, and how these states could be manip¬ 
ulated for practical applications of the TCIs. 

Despit e s ome strain-related studies in the 
literatur e ! 25 1 26 1 , we are not aware of a systematic in¬ 
vestigation of how spin-textures and other properties of 
Dirac states in the TCIs evolve under various symmetry 
breaking strains. So motivated, here we examine strain 
effects on the (001) surface of the Pbi_ x Sn x (Te, Se) 
TCIs. Our analysis is based on an effective four-band 
model, which is shown to capture all essential features 
of the topological surface states. The form of strain 
related perturbations on the Hamiltonian is clarified 
through general symmetry considerations, allowing us to 
delineate how characteristic features of the topological 
states in TCIs evolve under strains. 

Significantly, our analysis shows that certain mirror- 
symmetry-breaking strains induce hedgehog spin tex¬ 
ture with out-of-the-plane spin-polarization at the Dirac 
point. Such a spin texture of Dirac states has been re¬ 
ported previously in a 3D TI, but only in the presence of 
an exchange field 2 '. Our findings thus identify a possible 
new pathway for realizing TCI-based spintronics devices 
without requiring a TRS breaking fielcpS. Moreover, we 
show that the interplay of strains and applied Zeeman 
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field can be used to tune the Chern number of the surface 
ground state, which is a topological invariant character¬ 
izing a quantum anomalous Hall (QAH) or a quantum 
Hall insulator. 

Electron-electron interactions effects are ubiquitous 
in condensed matter systems, especially when charge 
screening is relatively poor as is the case generally for 
Dirac-like states in graphen d 29 * 30 ^ or TIs and TCIsPS In 
this connection, we discuss possible symmetry breaking 
orders generated under short-range repulsion U as a first 
step toward understanding correlation driven effects on 
Dirac states in the presence of applied strains. We nu¬ 
merically obtain the zero-temperature phase diagram as 
a function of U and the particle density n and delin¬ 
eate its evolution with increasing strength of strain. Our 
results not only give insight into the nature of compet¬ 
ing orders, but also build the foundation for what may 
be called “straintronics” applications driven by electron 
correlations. 

The paper is organized as follows. In Sec. II, we briefly 
discuss the effective four-band model for the (001) sur¬ 
face states in the low-energy regime, and its extensions 
for addressing strain effects in TCIs. Sec. Ill turns to 
consider strain effects on properties of the Dirac states, 
including the interplay between applied strains and per¬ 
pendicular Zeeman fields. These results then allow us to 
investigate in Section IV the possible interaction-driven 
competing orders in Dirac states in the absence as well 
as presence of strains. Finally, Section V comments on 
potential applications and implications of our study, and 
concludes with a summary of our results. 


II. EFFECTIVE FOUR-BAND k p MODEL FOR 
THE SURFACE STATES 

We start by reviewing the four-band k ■ p model for 
the (001) surface states in a TCI developed previously 
in Ref.1221 and discuss its generalization to account for 
effects of strains. On the (001) surface, the low-energy 
surface states can be viewed as two sets of interacting 
coaxial Dirac cones, originating from the interface be¬ 
tween the inverted bands of the TCI and the vacuum: 
one set is centered at X and the other at Y in the sur¬ 
face Brillouin zone (SBZ), see Fig. [TJ (b). Defining [110], 
[110], and [001] as x, y, and z directions [see Fig. [l] (a)], 
respectively, the effective model ffy around X must obey 
the following three essential symmetries which leave X 
invariant: the mirror reflection about the a,’,z-plane, the 
mirror reflection about the y^-plane, and time-reversal 
symmetry with the corresponding symmetry operations, 
i.e. 


M XZ H X (k X , ky)M xz - H X (k X , ky) , (1) 

M yz H x (k x ,k y )M y z = H x (—k x ,ky), ( 2 ) 

TH x (k x , fc !/ )T _1 = H x (-k x ,-k y ), (3) 

where k = ( k x ,k y ) is measured from X. 
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FIG. 1: (Color online) (a) Top view of the crystal structure 
of SnTe. Two types of atoms are drawn in different colors. 
x ( y ) denotes the direction along [110] ([110]). (b) Surface 
Brillouin zone of SnTe. Without strain, there are four gapless 
Dirac cones, protected by the mirror symmetry either with 
respect to FA' or FV. 


Informed by the orbital characters of surface states 
as revealed by first-principles calculations,^! one may 
choose, for instance in the case of SnTe, ||p 2 ,t;Sn), 
|p 2 ,F;Sn), |p 3 ;,t;Te), Ip^^Te)} as the basis states of 
H x because one of the interacting Dirac cones is p 2 - 
orbital (Sn) rich, while the other is p^-orbital (Te) rich. 
As a result, the symmetry operations can be represented 
by the matrices, M xz = —iT, 02 , M yz = — 1 X 31 , and 
T = —IT.Q 2 K, where Y, a p = a a G> sp with the Pauli 
matrices <x and s acting on orbital and spin spaces, re¬ 
spectively, and K denotes complex conjugation. These 
considerations restrict H x to the form 

H x = TOE30 + m'Yi 22 + k x (v\ x YiQ 2 + ^2x^20 + ^3x^32) 
+ ky( — VlyTlQl + V2y^l3 ~ V3yY,3l), ( 4 ) 

up to first order in |k|. 

For m '=0 and V2 X =V3 X =V2 y =V3 y =0, H x becomes 
block diagonal, reflecting the structure of the two un¬ 
derlying parent Dirac cones with two associated Dirac 
points, Ed\ = —m (p 2 -rich) and £^2 = to (p^-rich) at 
X [see Fig. [ 2 ] (a)]. A non-vanishing value of any of the 
parameters to ', V2 X , i’3x, V2 V and V3 V induces interaction 
between the two cones, leading to two immediate con¬ 
sequences: (1) The original Dirac points are shifted to 
£(0(2) = =F \/m 2 + m' 2 ] and (2) All degeneracies along 
the intersection of the two cones are lifted except for the 
two points of a time-reversal symmetric pair on the FX 
line, indicating the emergence of two child Dirac cones 
protected by the mirror reflection about the zz-plane [see 
Fig. [l] (b)]; hereafter, referred to as the low-energy Dirac 
points associated with £<^= 0 . 

Our effective four-band k ■ p model correctly capture 
two key features of surface states of TCIs. Firstly, 
for E > 0 topology of the constant energy contours 
changes from two separate Fermi circles at low energies 
to two concentric ellipses at high energies. This Lifshitz 
transition 15 indicates the presence of van Hove singulari¬ 
ties (VHSs) in the underlying electronic spectrum, which 
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FIG. 2: (Color online) (a) Band dispersion of the (001) surface of SnTe along high symmetry lines obtained by setting model 
parameters in Eq. as: a = 6.327 A, m = —0.3 eV, m! = —0.15 eV, vi x = V 2 V = V 3 X = V 3 y = 0 eVA, v\ x = —2.3 eVA, and 
Vi y = —6.5 eVA. Ed and Ed 1 ( 2 ) denote energies of the low-energy Dirac points and the two associated parent Dirac points, 
respectively. E VH s+<-) are energies of the van Hove singularities, (b) and (c) separately show spin textures of valence and 
conduction bands. Arrows and their lengths indicate directions and magnitudes of the in-plane spin polarization; the color map 
indicates the distribution of the out-of-the-plane spin polarization. Solid lines mark contours of constant energy. 


are located along the XM line at 

Ul 

E V hs+ = \m'\, k=(0,± —), (5) 

v ly 

if V 2 x, V 3 x, V 2 y , and v$ y are neglected. A similar situation 
arises for E < 0 with Evhs- = —Evhs+ • In all cases, 
we have a logarithmically diverging density of states pro¬ 
portional to In |)j| with u> and A representing an energy 
scale away from the VHS and a cutoff energy for sur¬ 
face states, respectively. Secondly, there is no out-of-the- 
plane spin component [Fig. [ 2 ] (b)-(c)] by using symmetry 
arguments in that the net out-of-the-plane spin polariza¬ 
tion, (£03) = (M xz M yz TT,o3T~ 1 M~^Mxz) = —(£03); 
yielding (£ 03 ) = 0. 

We turn next to generalize our four-band effective 
model to include low-energy states throughout the SBZ 
in the presence of inter-cone interactions. This can be 
done by noticing that the surface states in the vicinity of 
X and Y are related by C4 symmetry, so that the effec¬ 
tive model around Y can be obtained explicitly by a C 4 
rotation, 

Hy(k x , k y ) = CiHx{k y , — k x )C 4 1 , (6) 

where C4 = ao <8> e~ l T S3 and the p x orbital for Te atoms 
in the basis states is now replaced by p y orbital in H y. 
The total surface Hamiltonian then is 

■ff(ooi) = Ex (k x , k y ) ® Hy(k x , k y ). (7) 

Although F(ooi) describes an even number of Dirac cones 
at low energies like the case of a weak tP 2H311, note that 
the four Dirac points in a TCI do not locate at time- 
reversal invariant momenta. This key difference can lead 
to a rather different phase diagram in the presence of the 
electron-electron interactions 37 in a TCI, as discussed in 
Sec. IV below. 


III. STRAIN EFFECTS 
A. Gap opening and spin textures 

We now discuss effects on the surface states due to 
strains, which could be either intrinsic or extrinsic. Since 
the gapless surface states in the Pbi_ x Sn x (Te,Se) TCIs 
are mainly protected by mirror symmetries, a perturba¬ 
tion, which breaks one of these symmetries can be ex¬ 
pected to at least partly gap the surface spectrum, and 
modify the spin textures around the associated (massive) 
Dirac points. Although this is generally true, we will see 
below that this is not always the case due to other sym¬ 
metry considerations. For this purpose, we will carry out 
a systematic analysis of strain induced effects along the 
lines of Sec. II above. 

A general strain can be described in terms of a sym¬ 
metric strain tensor e^- = djUi, where u denotes the 
displacement held and i,j = x,y,z are chosen to coin¬ 
cide with the principal crystal axes as defined in Sec. II 
above. Any effective, strain-induced perturbation term 
corresponding to will transform in the same way as 
the strain tensor under time-reversal and mirror symme¬ 
tries. By examining the transformation properties of Y, a p 
up to first order in |k|, the possible resulting perturba¬ 
tions around X are listed in Table 0 

As expected, the strain induced perturbations in the 
first four columns of Table I can open gaps at both 
low-energy Dirac points along f A due to M xz broken 
symmetry. More significantly, the resulting spin tex¬ 
ture becomes hedgehog-like at low energies, similar to 
the spin reorientation phenomenon found on the surfaces 
of manganese-doped Bi 2 Se 3 thin Hlm^ 7 . 

As a concrete example, let us add the perturbation 
A £23 in Eq. <§> where A denotes the electron-phonon 
coupling strength. Clearly, in Fig. § a ), a gap opens at 
the Dirac point with its magnitude proportional to |A|. 
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FIG. 3: (Color online) (a) Band dispersion of (001) surface of SnTe along high symmetry lines with the strain perturbation 
AE23; parameters are the same as those used in Fig. [ 2 ] and A = 0.01 eV. (b) and (c) are spin textures of the conduction band 
in the presence of the strain and a perpendicular Zeeman field [modeled by adding H' = h z T ,03 in Eq. Q with h~ = 0.01 eV], 
respectively. Arrows and their lengths indicate the direction and the magnitude of in-plane spin polarizations; color maps show 
the distribution of the out-of-the-plane spin polarization. Solid lines mark constant energy contours. 


TABLE I: Effects of perturbations induced by applied strains. The first row shows the form of each perturbation in terms of 
a 4 x 4 matrix in the basis of the effective model around A'. A and v are the coupling strengths. The next two rows indicate 
if a perturbation breaks certain mirror symmetry. The fourth and fifth rows represent the possible consequences. The last row 
shows the corresponding strain tensor with which each perturbation can couple. 


Perturbation 

AE23 

vk x E03 

vk T An 

vk x E 33 

AE10 

vkx E12 

AE21 

vkx^oi 

vkx E13 

vk x E31 

Break a;z-mirror symmetry 

/ 

/ 

/ 

/ 

X 

X 

/ 

/ 

/ 

/ 

Break j/z-mirror symmetry 

X 

X 

X 

X 

/ 

/ 

/ 

/ 

/ 

/ 

Open gaps at Dirac points 

/ 

/ 

/ 

/ 

X 

X 

X 

X 

X 

X 

Induce out-of-plane spins 

/ 

/ 

/ 

/ 

X 

X 

X 

X 

X 

X 

Coupled strain tensor 

Syz 

£yz 

£yz 

Syz 

£xz 

£xz 

S X y 

S X y 

S X y 

S X y 


The spin-texture of the conduction band of the massive 
Dirac cone is also depicted in Fig. [3])b). As the absolute 
value of the eigenenergy |E| becomes smaller (be., closer 
to the bottom of the upper cone), the induced out-of-the- 
plane spin component, 

_ 2m’Av lx k x 

[ 03) ~ E[E 2 - (m 2 + m ' 2 + v\ x k 2 x + v 2 ly k 2 y + A 2 )] ’ 

( 8 ) 

becomes larger, where we have set v^ x = V 2 y = v^ x = 
Vs y = 0 for simplicity. Several points should be noted 
here as follows. Firstly, since the strain perturbation 
does not break TRS, the spin texture for the two mas¬ 
sive cones around X must form time-reversed partners. 
This can be inferred from Eq. ([ 8 ]) by changing the sign of 
k x . Secondly, the out-of-the-plane spin component of the 
lower cone is basically opposite to that of the upper cone 
at the same |k|. Finally, the whole spin texture resulting 
from the strain differs from the case where the hedgehog¬ 
like texture is induced by the perpendicular Zeeman field 
when both massive Dirac cones around X are considered, 
see Fig. [3])c) . 

Other strain perturbations in Table [IJ which break ei¬ 
ther the yz-mirror symmetry or both the mirror symme¬ 
tries, neither open a gap at the Dirac point nor induce 
an out-of-the-plane spin component. The robustness of 


this gapless Dirac point originates from an underlying 
symmetry: the former type of perturbation is due to the 
presence of M xz symmetry, while the latter type (shear 
deformation) is protected by a “space-time” symmetry, 
Ci t — *E 31 A", with a rotation C 2 followed by a time- 
reversal operation. The space-time protection now al¬ 
lows persistence of gapless Dirac points sitting at generic 
k points away from the mirror line TA. This special fea¬ 
ture has been noted previously in Refs l 2 ol 2 (H and it leads 
to a new type of TCpO. 

Vanishing out-of-the-plane spin polarization can be 
proven via symmetry arguments. For instance, consider 
the surface states around X under the strain perturba¬ 
tion AE 10 . The corresponding Hamiltonian now reads: 
Hx = Hx V ASio- Assuming 141) is an eigenstate of L/y, 
the xz-mirror symmetry guarantees that 

0 = |[Eo2,£jf]|*> = 2*u ll ,fc I ,<tf |E 03 |tf). (9) 

Thus, the out-of-the-plane spin component ('F|E 03 |4/) 
vanishes everywhere around X. Other perturbations can 
be analyzed in a similar manner. 

There other cases deserve some comment. The first 
two cases involve a uniform expansion and a uniaxial 
stretch (C 4 breaking), which correspond to e xx + e yy 
and e xx — Eyy , respectively-®. Because both these cases 
respect TRS as well as the mirror symmetries, their 
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net effect is only to renormalize the parameters in the 
original H x - Consequently, these perturbations do not 
open a gap at Dirac points, although positions of the 
Dirac points could shift in opposite directions along the 
TA line. The third case is an experimentally observed 
ferroelectric-like distortion 211 , in which two kinds of atoms 
are displaced along a certain direction in an opposite 
manner. Denoting the displacement vector d = (d x ,d y ), 
the non-vanishing component of this distortion, d x ( d y ), 
preserves the mirror symmetry with respect to the prin¬ 
cipal axis x ( y ), but breaks the rotation and the mirror 
symmetries along a perpendicular direction. To the ze¬ 
roth order in k around A, the perturbation due to such 
distortion can be straightforwardly shown to have the 
form: 


strain and Zeeman fields. An applied Zeeman field, which 
respects C 2 rotational symmetry around the z-axis, opens 
gaps with the same sign in the pair of Dirac cones around 
A, while an applied strain held, which breaks C 2 rota¬ 
tional symmetry, would induce gaps with opposite sign. 
This observation provides the foundation for tenability 
of the Chern number via the interplay of applied strain 
and Zeeman fields, suggested in Refill using a simplified 
two-band model for each Dirac point. 

As a concrete example based on our four-band model, 
consider an applied Zeeman held along the z direction 
along with a strain held which breaks M yz but preserves 
M xz around Y on the (001) surface. From symmetry 
considerations, the strain and Zeeman helds around X 
and Y become coupled as 


H' f = Xp x d x YiiQ + Aj7 y dyE23, (10) 

where Ap x and Ap y denote the coupling strengths. In 
fact, the former term is similar to the effect of strain e xz , 
preserving the gapless Dirac points, while the latter term 
is similar to the effect of strain e y2 , resulting in opening 
a gap instead. 



H x + h z Y, 03 , 

Hy + h z E 0 3+A s E 23 , (11) 

where h z (V),) is the held strength of Zeeman (strain) 
held. Note that for a given h z , all four low-energy Dirac 
points open gaps in the absence of an applied strain [see 
Fig. Ftla)]. When a non-vanishing V s <C \h z is intro¬ 
duced) the gaps of the two massive Dirac cones around 
Y evolve in an opposite manner in the sense that one 
cone increases while the other decreases. This is con¬ 
sistent with a picture in which the strain held induces 
out-of-the-plane spin components around the two mas¬ 
sive Dirac cones which form time-reversed partners with 
respect to their spin textures [see Fig. |3][b)], and there¬ 
fore, respond oppositely to the existing h z held. With 
increasing strength of V s , the gap in the Dirac cone with 
a decreasing gap continues to decrease further, becomes 
gapless at V s = \h z , and reopens inverted again, as 
shown in Figs. |4|b)-(d). This band inversion indicates 
that the system undergoes a topological phase transition. 

The topological nature of a ground state with bro¬ 
ken time-reversal symmetry can be characterized by the 
Chern number. A non-zero value of the Chern number, 
C, given by 


FIG. 4: (Color online) Evolution of band gaps in the Dirac 
cones under strain perturbation I 4 E 23 with a fixed Zeeman 
field, where h z = 0.1. The strength of the strain perturbation 
is given by V s = (a) 0, (b) 0.03, (c) 0.05, and (d) 0.08. From 
(b) to (d), one of the Dirac cones around Y closes and reopens 
again due to a band inversion. The model parameters are the 
same as those used in Fig. [2] 


B. Tuning the Chern number via the interplay of 
strain and a Zeeman fields 

Like a perpendicular Zeeman field, the strain field on 
the surface of a TCI may not only give rise to spin reori¬ 
entation, but it may also induce a charge gap in the Dirac 
cones. However, there are key differences in the effects of 


C = 


-1 
7r 


E E 

nGocc. m^n 


Im[(z 


dk x \ u m)(Um\Qk y l Ura )] 

{E n - E m y 


( 12 ) 


where the summation over momentum space is implicit, 

2 

indicates a finite Hall conductance a = C^ 1 which can 
be obtained by integrating the Berry curvature of the 
wave functions of the occupied bands over the momentum 
spac^Sl. Fig. [ 5 ] presents the Chern number of the system 
as a function of the Zeeman held strength h z and the 
strain strength V s within the framework of our four-band 
model. In the presence of a Zeeman held, by varying the 
strength of strain from small to large values is seen to 
drive the system from C = 2 to C = 1, demonstrating 
the tunability of the system to the strain held. Note 
that in our consideration here we have implicitly assumed 
that the sample is thick enough so that the hybridization 
between the top and bottom surfaces of the sample can be 
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neglected; with h z ^ 0 the bottom surface will contribute 
another Chern number C = 2, which is not shown in 
Fig-Ill for simplicity. 



FIG. 5: (Color online) Chern number for our four-band model 
of a TCI as a function of 14 and h z . Model parameters are 
fixed as: m=- 2, m'=-l, vi x =-l, vi y =-1.5, with other parame¬ 
ters taken to be zero. The phase boundary (blue line) is given 
by h z = 214- 


IV. EFFECTS OF ELECTRON-ELECTRON 
INTERACTION 


In discussing electron-electron interactions, we con¬ 
sider short-range repulsive interactions between electrons 
on the (001) surface of a TCI, and focus on delineat¬ 
ing how the phase diagram evolves in the presence of a 
strain, which breaks one of the mirror symmetries, say, 
M yz , based on the the interacting Hamiltonian: 

H = Y,ft(k)H {001) iP(k) (13) 

k 


where 


^ f (k) = [ C t 

j 


x,p z t( k )’ c jc, P JV’ c kp*tM’ c kpJ k) ’ 


r, Pz t 


(k), 


r,pA 


( k )’ c L„ t (k) ’ 


4, P ^( k )K 14 ) 


with the subscript X (Y) denoting the momentum point 
involved in the k expansion. For H( 0 oi), we have added 
the strain perturbation and hence, F/(ooi) = H y ® (Hy + 
V s E 2 3 ). The interaction term, H int = Hu + Hy-, Hjj 
and Hy denote repulsive contact interactions between 
like and unlike orbitals, respectively, with 

Hu = U d 2 r Y ^t( r )^t( r )> 

V=Pz,Px,Py 

H V = j f d 2 r ^n^(r)n v (r), (15) 

1 


where h v ( r) = E>^( r ) = Eq, s (Q = 


X,Y, r) = p x ,p y ,p z ) [c Qj7)s (r) are defined in Eq. (|16j) 
below]. 

Since we are treating only the low-energy, long- 
wavelength physics, it is reasonable to consider field op¬ 
erators CQ r/a (r), which vary slowly on the scale of the 
lattice constant. This can be done by Fourier transform¬ 
ing the operators: 


cq.^r) = Y e lK r c Q ,, s (K) = e lQ r c Q , I)s (r), (16) 

K 


where K = Q + k, expanding from the origin (0,0). In 
terms of these field and density operators, and the iden¬ 
tities given in the Appendix, the full Hamiltonian H can 
be rewritten as 


H = / d 2 r&H m) * 


■H, 


(17) 


where (ET (r) has the same form as Eq. (14) except that 
the k-dependence is now replaced by r. 

In order to determine the ground state of the system 
as a function of the interaction strength and fermion den¬ 
sity, we employ the self-consistent mean-field (MF) ap¬ 
proach, which is expected to be reliable as long as the 
interaction strength is much smaller than the bulk band 
gap. Taking U = 14 as a representative case, we de¬ 
couple the Hi n t into bilinear fermion terms. After some 
straightforward but tedious algebra, we obtain 52 order 
parameters as well as 12 renormalized band parameters 
as detailed in the Appendix, which can be generally ex¬ 
pressed as 


O ia p = ('F t E 7a/3 'F) (18) 


where a summation over momentum and occupied states 
is implicit, and = r 7 0 £ a/ g, with the Pauli matri¬ 

ces, t 7 , acting on XY pseudo-spin (valley) space. Note 
that, here we consider only the symmetry breaking or¬ 
ders in the particle-hole channel: superconductivity in 
the particle-particle channel due to weak onsite repulsion 
can only be achieved via beyond MF treatment, usually 
with exponentially small transition temperature 41 . 

In this connection, we define the “particle density” at a 
given chemical potential p as the density deviation from 
the total particle density where p = E^ = 0, namely, 


n = 


rE( k)= M 
/ k<kA 


d 2 k 

(2^F 


l*Eo(k)=0 
' k<kA 


d 2 k 

W’ 


(19) 


where i?o(k) (E( k)) is the eigenenergy of the non¬ 
interacting (MF-decoupled interacting) system; k a is a 
cutoff momentum, which is chosen such that our results 
are insensitive to its magnitude. Note that by any in¬ 
duced valley “polarization”, O 300 , we mean 

r E o(k)=0 ^ 2 ^, 

O 300 = (d'^E 3 oo4 / ) — / . 2 (E 30 o). (20) 

J k<kA \^J 
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FIG. 6: (Color online) U versus n phase diagram without 
strain. nvHS denotes the particle density at Evhs+- Solid 
and dashed lines mark phase boundaries of continuous and 
abrupt phase transitions, respectively. 


That is, the “valley polarization” due to external strain 
in the absence of interactions is subtracted in defining 

0300- 

For a given U and n, the ground state can then be 
determined by minimizing the MF free energy with re¬ 
spect to various order parameters using the form detailed 
in the Appendix. We numerically solve the resulting set 
of coupled equations self-consistently to obtain the zero- 
temperature U-n phase diagrams both with and with¬ 
out the applied strain. Note that, in order to make 
the underlying physics more transparent, we use the fol¬ 
lowing model parameters in computations: m = —2, 
m! = —1, v\ x = —1, vi v = —1.5, and a cutoff energy, 
A = Ef)(k\) = 2 \Jm 1 + to' 2 , which determines k\. How¬ 
ever, we expect our results to be generic and relevant 
more generally, detailed effects of material-specific pa¬ 
rameters notwithstanding. 


A. System without strain 


At n = 0, the “p z -CDW” phase appears when the 
interaction strength U is larger than the critical value 
U c ~ 2.7. This phase involves a finite value of Ol °o+°i3o ; 
it is associated with broken translation symmetry, while 
the TRS and both the mirror symmetries remain intact 
[see schematic Fig. [7))b)], so that each Dirac cone remains 
gapless. The existence of a critical value of U can be un¬ 
derstood by noticing that in the non-interacting limit, 
there are only four Dirac points with zero density of 
states at the Fermi level. The necessity of a threshold 
value of U has also been predicted theoretically in 2D 
systems with linear energy dispersion, such as graphene 
and the surfaces of 3D weak TltPHHl. 

Upon electron doping, the system assumes either the 
FM phase [see Fig.jTja)] at relatively high particle densi¬ 
ties, or the p z -CDW phase [see Fig. [7j[b)] at low particle 
densities. The FM phase is associated with broken TRS 
as well as both the mirror symmetries simultaneously, 
resulting in gaps to open up at all the Dirac points, al¬ 
though the corresponding ground state is still metallic 
with finite p. Note that the critical U, above which the 
system enters into the gapless p z -CDW phase, is less than 
the critical value for the FM phase, indicating a favorable 
free energy gain compared to the competing FM phase 
at low doping. 



y | - 

|U0| 

U . 


tie i 

r» ® I 




FIG. 7: (Color online) Schematic plots illustrating (a) the 
ferromagnetic state (FM), and (b) the p z -orbital rich, (n, n) 
charge density wave state (p z -CDW). Different colors label 
the two distinct atoms involved; arrows depict the net out-of- 
the-plane spin polarization. 


We begin by considering a system without strain, i.e., 
V s = 0. In the order parameter space we explored in this 
case, the leading orders (aside from the metal phase) are 
basically the ferromagnetic state (FM), O 003 , and the 
p z -orbital rich (n,n) charge density wave (p z -CDW): 

Ql0 °f 130 = (^(- 1 )^[n Pz ,t(r) + n Pa4 (r)]), (21) 

r 

where the length scale is in units of the lattice con¬ 
stant. The zero-temperature U-n phase diagram is shown 
in Fig. [6] The p z -CDW phase corresponds to a non¬ 
vanishing Ol oo+ Ol 3o ^he i ea ding piece), and the FM 
phase represents a non-vanishing value of O 003 (the lead¬ 
ing piece). Finally, in the metal phase, all symmetries 
are preserved with vanishing values of all 52 order pa¬ 
rameters. 


We further note several points in connection with the 
phase diagram of Fig. [6] as follows: (1) The spin density 
wave (SDW) phase, as anticipated usually in Dirac sys¬ 
tems like graphene, does not occur here at n = 0 for large 
U. This could be attributed to the non-trivial orbital 
textures and strong spin-orbit coupling in our case; (2) 
We find that the dynamically generated strain-like order, 
which could gap the spectrum, is relatively disfavored as 
the system tries to preserve the combined C 2 T symmetry, 
although such a phase becomes favored in the presence of 
external strain; (3) If we suppress inter-valley scattering, 
our preliminary results suggest that the CDW phase is 
also suppressed^! and finally, (4) As we pointed out al¬ 
ready in Sec. II above, the density of states at an energy 
scale u! near a van Hove singularity diverges like In |w|, in¬ 
dicating propensity for the occurrence of more symmetry 
















FIG. 8: (Color online) U versus n phase diagram in the pres¬ 
ence of a moderate strain. Lines mark boundaries between 
various continuous phases. 


breaking orders (including those in the particle-particle 
channel). The present mean-field treatment would then 
be inappropriate, even in the weak-coupling limit. 


B. System under moderate strain 

Here we consider an applied strain with a moder¬ 
ate coupling strength, V s = \m'\ [of order 0{Evhs+)\, 
which breaks the yz -mirror symmetry. The resulting 
zero-temperature phase diagram is shown in Fig. [8j In 
contrast to the case without strain, a coexisting phase ly¬ 
ing between the p z -CDW and FM phases is now seen to 
emerge. Suppression of the p 2 -CDW phase with strain 
can be anticipated on physical grounds. The reason is 
that any (7r,7r) order requires good nesting between the 
valleys, but the strain weakens the nesting tendency by 
opening gaps at Dirac points along the yz-mirror line. 
Notably, in the strained system, spin polarization in the 
FM phase is mainly contributed by the p x orbital due to 
the broken mirror symmetry, which lifts the p y orbital to 
higher energies. 


C. System with large strain 

When the coupling strength of the applied strain V s 
becomes the largest energy scale in the system, all other 
degrees of freedom can be ignored at the Y valley as 
states around Y are all gapped out, leading to the phase 
diagram of Fig. 9} All (n, 7r) orders are seen to disap¬ 
pear, and only the FM phase survives for U > U c (n) at 
a given density n. The presence of the concave boundary 
region separating the metallic and FM phases in Fig. 8 
can be understood from energetics: At any given den¬ 
sity in this region, the chemical potential lies within the 
gap induced by the FM order at 2 , yielding a gain in 
free energy and thus lowers U c . As the particle density 


FIG. 9: (Color online) U versus n phase diagram in the case 
of large strain. Solid line marks the boundary between the 
two continuous phases. 


increases, the chemical potential is eventually unable to 
locate within the gap and this argument fails. Interest¬ 
ingly, at n = 0, although the ground state is insulating 
with a small gap due to TRS breaking, the system does 
not exhibit quantum anomalous Hall effect. This should 
be contrasted sharply with the case of a Zeeman field ap¬ 
plied perpendicular to the (001) surface of a TCI in the 
non-interacting limit, where such a broken TRS leads 
to non-vanishing Chern number as shown in Sec. IIIB 
above. 


V. DISCUSSION AND CONCLUSION 

Some of the strains discussed in this study could 
be realized experimentally through chemical doping. 
Since PbSe and SnSe assume cubic and orthorhom¬ 
bic structured, respectively, strains or lattice distor¬ 
tions could likely be manipulated by adjusting chemi¬ 
cal compositions. In fact, a ferroelectric-like lattice dis¬ 
tortion has been observed in the TCIs via STM/STS 
spectroscopy 21 . Mechanical tuning could provide another 
approach. For example, a piezoelectric layer could be 
grown on the surface of a TCI to allow manipulation of 
strains via electric field as proposed in Refl25) 

As to practical applications, we have shown clearly 
that the Chern number of the surface states can be tuned 
through the interplay of applied strain and Zeeman fields. 
This provides a potential new handle for controlling 
the topological conducting edge channels for low power¬ 
consuming, next generation electronic technologies. Fur¬ 
thermore, our analysis shows that in the presence of 
electron-electron interactions, ferromagnetism could be 
turned on/off via strain fields, which, offering a novel 
pathway toward spintronics applications. 

In summary, we have systematically investigated the 
effects of strains, which break either one or both mir¬ 
ror symmetries without breaking time-reversal symme- 
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try, on the (001) surface of the TCIs. Under mirror- 
symmetry breaking strains, the low energy Dirac cones 
not only become massive, but also develop hedgehog¬ 
like spin textures. We show that the Chern number of 
the surface states can be tuned via applied strain Helds 
in the presence of a Zeeman field. Finally, we delin¬ 
eate the competing orders that can result from varying 
strengths of short-range (repulsive) electron-electron in¬ 
teractions at the mean-field level, and show that the re¬ 
sulting correlated phases are amenable to control through 
strain fields. The tunability of the interacting as well as 
the non-interacting electronic structures and topological 
states of the TCIs via strain fields and their interplay 
with Zeeman fields as revealed by our study suggests 
new pathways for developing spintronics and other ap¬ 
plications platforms based on the TCIs. 
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Appendix A: Derivation of mean-field theory in Sec. 

IV 


Here we present further details of the MF theory used 
in Sec. IV. In terms of the slowly varying held operators 
defined in Eq. (16), the relevant number operators are 
given by 


hu t( r ) — c j?,p*t c - x >*t 


c y, Pz ^ c y,pA 


+ (~l) x+y (c\ >Pz ^ c Y,p z \ + C Y,pA C X 
’WM = Jx iPA Cx,p z \, + Jy,p z i c y, P z I 

+ (~ l ) X+V ^X,p z \. c Y,p--l + c F,p,4. c X,pn)> 


Op 


u( r ) = 


W( r ) = 
W( r ) = 


c t 
x, Px t 

J 

x, Px i 


- J 


°Py 


4 .( r ) = 

s ( r ) = 
= ( r ) = 


( r )cA,p xt ( r )’ 

(r)cA,pa( r )> 

-Y, Py ^ r ) c Y,Py t( r )> 

c Y, Py i( r ) c Y,P y l( r ), 
npMr) + n Pzi {r) : 


% p* t 


(r) 




°Py 


( r ) = ™P H t( r ) +«P„ 4 .( r )- 


(Al) 


We illustrate our MF treatment with the example of 
J d 2 r Uh Px ^(r)n Pxj A r )- For this purpose, we decouple 
bilinear fermion operators composed of A or B in the 
particle-hole channel into A(B) + B(A) — (B)(A) + • • •, 
where • • ■ represents fluctuations away from the mean- 
held values, which are neglected. The condensation en¬ 
ergy, — (B) (A), for f d 2 r Un p f(r )n p ,|,(r) is thus given 
by 


+ S 300 + £303 ~ £330 — £ 333 )’!' 

g 

= ^(0000 + 0003 — 0030 ~ 0033 
+ 0300 + 0303 “ 0330 ~ 0333) 5 


(A3) 


where O ia p was dehned in Eq. (18) and S' is a unit 
area. Notice that O 0 00 is the particle density dehned 
in Eq. ©. Whe n the similar procedure is applied to 
other terms in Eq. (A 2 ), it turns out that —(B) (A) term 
becomes: 


^^64 t(O° 00 — ^ 030 ^ 300 ~ 033o ) 2 

— (0003 — 0033 + 0303 ~ 0333 ) 2 ] 
+ 5 '^ [(0001 — 0031 + 0301 — 0331 ) 2 

+ (0002 — 0032 + 0302 — 0332 ) 2 ]- 


(A4) 


This allows us to obtain the MF decoupled 
f d 2 T Un Px (r)n Po .4 (r) as 

J ^ 2r { 32 K^ 000 — ^ 03 ° + ^ 30 ° — ^ 33 °) 

x 4't(Eooo — £030 + £300 — Ssso)'! 3 
— (0003 — 0033 + 0303 — 0333) 
x 'h^(£o03 — £033 + £303 — Ssss)'! 3 ] 

32 ^(O° 01 — ^° 31 "b < ^ 301 — 033i) 


x 4't(E 0 oi — S 031 + S 30 i — E 331 )T 
+ (0002 — 0032 + 0302 — 0332 ) 

X 4lt(Eoo2 — £o32 + £302 — £332)'h] j 

~^64^0 OOCI ~ ^ 030 ^ 30 ° — ^ 33 °) 2 

— (0003 — 0033 + 0303 - 0333 ) 2 ] 

+<5 [(0001 — 0031 + 0301 — 0331 ) 2 

+ (0002 — 0032 + 0302 — 0332) 2 ]- (A5) 


By applying the same trick for the other interaction terms 
in Eq. (15), we finally find 64 MF undetermined param¬ 
eters (summarized in Table. 0 , and obtain the MF free 
energy per unit area at zero temperature as: 


“ U (4c, Px t c x,p^)( c X, Pol i Cx ’P^) 

+ U ( 4 c, Px t C X’P*i){ C X, Px l C X,p*t)’ ( A2 ) 

where the r-dependence is implicit. In fact, 
( C A p f c A,p a .t) can b* e further rewritten as 


S 

8 


^ t (£ 0 oo + £003 


— £q30 


— £033 


F= lj J d 2 r¥ (if(ooi) + H% t F )$ ~ny + E con , (A 6 ) 

where is the MF decoupled interaction Hamilto¬ 

nian for Hi n t and E con is the MF condensation energy. 
The ground state can then be determined by solving the 
set of coupled mean-held equations resulting from mini¬ 
mizing F with respect to all possible order parameters. 
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TABLE II: Summary of MF parameters 


Type 

MF parameters 

Renormalized 

O 030, 

C?022, 

a 0322 ^ 0021^ 

0321, 6 

band parameters 

Ooio, c 
O 330/ 

O310 

Oooo 

, C Oo23, C 0323, C 

O300/ 


O002, 

0032, 

O001, 

O 031, 

Oo03, 

O 033, 

FM 

O013, 

O313, 

O 020, 

O 320, 

Ooil, 

O311, 

O012, 

0312, 

O 303, 

O333, 

O30I, 

O331, 


O 302, 

0332 






O100, 

Ol30, 

O110, 

Ol23, 

Ol22, 

Ol21, 

CDW 

O220, 

0213, 

O212, 

O2II, 

O 203, 

0233, 


O 202, 

0232, O 201, 0231 




Ol03, 

Ol33, 

0102, 

O132, 

0101, 

O131, 

SDW 

O113, 

Ol20, 

Oin, 

O112, 

O 223, 

O 210, 


C?221, 

0222, O 200 , 0230 




a In the absence of strain, if O022 7^ O322 , this is an order param¬ 
eter for spontaneous mirror symmetry breaking. 

b ln the absence of strain, if O021 7^ — O321, this is an order pa¬ 
rameter for spontaneous mirror symmetry breaking. 

c In the absence of strain, if one of O010? O310, O023 and O323 is 
nonzero, this would be an order parameter for spontaneous mirror 
symmetry breaking. 

rf In the absence of strain, if either O300 or O330 is nonzero, this 
would be an order parameter for spontaneous C4 symmetry break¬ 
ing. 



